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ABSTRACT 

We show that the spin period of the white dwarf in the magnetic CV EX Hydrae 
represents an equilibrium state in which the corotation radius is comparable with 
the distance from the white dwarf to the inner Lagrange point. We also show that a 
continuum of spin equilibria exists at which Pgpin is significantly longer than ~ O.lPorb- 
Most systems occupying these equilibrium states should have orbital periods below 
the CV period gap, as observed. 

Key Words: accretion, accretion discs - binaries: close - stars: individual: EX Hydrae 
- stars: magnetic fields. 



1 INTRODUCTION 



Cataclysmic variables (CVs) are semi-detached binaries in 
which a low-mass, main-sequence secondary star transfers 
msiss to a white dwarf. In a significant number of systems 
the white dwarf has a magnetic moment strong enough to 
channel the accretion flow on to restricted regions of its sur- 
face; such systems are known as magnetic CVs. Observation- 
ally these systems are characterized by the presence of two 
coherent periods Porb, -Pspin, corresponding to the orbital ro- 
tation of the binary and the white dwarf spin respectively. 
Sometimes the beat period 



(P~ 

\ spin 



(1) 



is seen in addition to or even in place of Pspin. Magnetic 
CVs are interesting because the white dwarf spin rates give 
insight into the angular momentum flows within the binary. 
These flows are more complex than in the analogous systems 
(pulsing X-ray binaries) in which the accretor is a magnetic 
neutron star, simply because the white dwarf magnetic mo- 
ments fj-i 10^^ G cm^) are much larger than those of the 
neutron stars (^i < 10^° G cm"^). Accordingly the magne- 
tospheric radius in a magnetic CV is a significant fraction 
of the binary separation a, whereas it is far smaller in the 
neutron-star case. As an illustration of this, the white dwarf 
spin can be locked to the binary rotation (Pspin ~ Porb) in 
some magnetic CVs (called AM Her stars), while Pspin is 
always considerably smaller than Porb in neutron-star sys- 
tems. 

Historically the neutron-star systems were discovered 
first, and simple theories of the angular momentum flows 
were developed for them (e.g. Ghosh & Lamb, 1979). There 
has been a natural tendency to transfer these theories bod- 
ily to the white dwarf case, even though the larger mag- 
netic moments make this very questionable. In particular, 
neutron-star magnetospheres are so small that Roche lobe 
overflow must always lead to the formation of a Keplerian 
accretion disc in the standard way (e.g. Frank et al., 1992). 
By contrast, a white dwarf magnetic field can have a strong 



infiuence on the accretion flow at all radii, making it non- 
Keplerian even though it may surround the white dwarf. 
This in turn tends to lead to an equilibrium spin rate deter- 
mined approximately by the condition (King, 1993; Wynn 
& King, 1995) 



Pco ^ Pcirc . 

Here 

/ GMlPgpin^ ^ 



R. 



\ 47r2 



(2) 



(3) 



is the corotation radius, at which matter in local Kepler ro- 
tation about a white dwarf of mass Mi corotates with the 
magnetic held lines, and Pcirc is the circularization radius, 
i.e. the radius at which the Kepler specific angular momen- 
tum equals that of matter accreting through the inner La- 
grange point Li. If the latter is a distance b from the white 
dwarf we can find an approximation for Pcirc by equating 
the specific angular momentum ~ 6^(27r/Porb) at Li to the 
specific Kepler angular momentum at Pcirc (this neglects 
the effects of the companion) 

1/3 



Rc 
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(4) 



where M — Mi + M2 is the total binary mass. Numerical 
fits to Roche geometry give 



- = 0.500 -0.227 log 

a Ml 



(5) 



(Warner 1976). Hence combining (^ - |^ gives a spin equi- 
librium 



p . 

J spin 
Porb 



M 

aTi 



0.07 



(6) 



for typical mass ratios (cf King, 1993, Wynn & King, 1995). 
The numerical simulations performed by Wynn & King 
(1995) used the magnetic drag prescription introduced by 
King (1993). Here the flow is assumed to be in the form 
of diamagnetic blobs, the local magnetic held B = ni/r^ 
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acting on them through a surface drag. This produces an 
acceleration 



-fc[v - Vf]x, 



(7) 



where v and Vf are the blob and field velocities (i.e. V{ = 
27rr/Pspin), the suffix 1. denotes the velocity components 
perpendicular to the field lines, and k 
drag time scale given by 



imagi with tmag the 



tn 



(8) 



(Drell, Foley and Ruderman, f965). Ifere pb, are the blob 
density and length scale, and ca is the Alfven speed in the 
inter-blob plasma. Thus as the blobs move through the mag- 
netosphere they exchange orbital energy and angular mo- 
mentum with the white dwarf via the drag term. The blob 
parameters pb , Ih are to some extent arbitrary, but limited 
for example by physical conditions near Li. Numerical sim- 
ulation of the flow including the full Roche potential rather 
than the simple approximation examined the regime in 
which 

tmag(-Rcirc) ^ ^dyn (-Rcirc) , (9) 

or alternatively 

tmag(I/l) > tdyn(-Ll), (10) 

where t^yn ~ {r^ /GMi)^''^ is the dynamical time scale. 
These simulations give similar but slightly larger values of 
fspin/^orb to the approximation 

Equation ^ reasonably describes the relation between 
Papin and Porb in most non-synchronous magnetic CVs ('in- 
termediate polars', or IPs). The corresponding accretion flow 
is sometimes referred to as 'discless'. This is a rather mis- 
leading name, suggesting for example that the flow cannot 
surround the white dwarf, although simulations show that 
this is possible (Fig. 1). Observations suggesting for exam- 
ple that matter is fed to the white dwarf at all spin phases 
are sometimes taken as evidence against this type of accre- 
tion. However the fundamental property of this flow is not 
the absence of matter surrounding the white dwarf, but sim- 
ply the non-Keplerian nature of the velocity field. This in 
turn reflects the pervasive influence of the magnetic field on 
the accretion dynamics. The divergence from Keplerian flow 
is a manifestation of the complex angular momentum flow 
within the system, where the angular momenta of the blobs 
may be transferred to the primary star via the field or accre- 
tion, recaptured by the secondary star or lost to the binary 
entirely. This contrasts with the case of a Keplerian accre- 
tion disc where the angular momentum transferred through 
Li is passed back to the binary orbit via tidal forces at the 
outer edge of the disc. 

Since the IPs cover a large range of magnetic field and 
orbital period there are several exceptions to the equilib- 
rium represented by One would expect systems with 
Keplerian discs to have Pco < Pcirc and thus a smaller ratio 
Pspin/Porb thau giveu by (^. This appears to be the case in 
at least two systems. In GK Per (Pspin — 381 s, Po^b = 48 hr) 
the orbital period is so long that the magnetospheric radius 
is small compared with Pcirc , so a Keplerian disc will form in 
the standard way. In DQ Her (Pspin = 71 s, Porb = 4.65 hr) 
the very short spin period presumably indicates a rather 
low magnetic field, leading to the same conclusion. How- 



ever, the still shorter spin period of AE Aqr (Pspin = 33 s, 
Porb ~ 9.88 hr) does not imply a Keplerian disc (Wynn, 
King & Horne, 1997, and references therein); here the mat- 
ter transferred from the secondary star is apparently being 
centrifugally ejected, showing that disc formation depends 
on Pspin as well as the magnetic moment fii (cf. equation 7). 

In addition to the systems with Pspin/Porb smaller than 
implied by (^ there are two IPs with a much larger ratio: 
EX Hya has Pspin = 67 min, Porb = 98 min, and RX1238-38 
has Pspin = 36 min, Porb ~ 90 min. This clearly does not 
fit the 'IP' spin equilibrium (^), still less the equilibrium 
expected for accretion from a Keplerian disc. A spin period 
far from equilibrium would imply that we are seeing these 
systems in a very short-lived phase: in the case of EX Hya 
the spin period derivative is Pspin = —3.8 x 10~^^ s s~^ (eg. 
Warner, 1995), with a spindown timescale of ~ 10® y, much 
shorter than the orbital evolution timescale (^10* y). It is 
far more likely that the spin is currently close to equilibrium, 
but that this equilibrium differs from that represented by 
(^) or the Keplerian case. In this paper we shall show that 
EX Hya represents a new type of spin equilibrium for IPs 
with short orbital periods. 



2 SPIN EQUILIBRIUM IN EX HYDRAE 

Most IPs so far discovered have orbital periods on the long 
side of the CV period gap 2 hr < Porb ^ 3 hr. EX Hya's short 
orbital period is therefore very suggestive: if its magnetic 
field is similar to systems conforming to ^ we would expect 
the smaller orbital separation and lower mass transfer rate to 
make the effective magnetospheric radius comparable with 
the distance to the secondary star. Here magnetic forces will 
dominate the motion of the gas blobs near Li and 



tmag(il) i^tdyn{Ll). 



(11) 



Setting Pb and lb to be equal to the expected stream density 
we find t 



10"^ g cm"^) and width (- 



mag(il) 



10 cm) at Li respectively, 
tdyn(ii) for Mi(crit) - 3 x 10^^ G cm^ 
(cf fig. 6). For magnetic moments larger than this we might 
then expect to see the equilibrium condition (^ (which leads 
to eq. p|) replaced by 



Pco ~ b, 

leading (from |^, |^ to 

3/2 



Porb 



0.500 - 0.227 log 



Ma 
Ml 



3/2 



(12) 



(13) 



For the mass ratio M2/M1 = 0.19 (Hellier, 1996) appropri- 
ate to EX Hya we find Pspin/Porb ~ 0.54, similar to the 
measured Pspin/Porb = 0.68. However this simple treatment 
completely neglects the full Roche potential. Since the angu- 
lar momentum exchange determining Pspin/Porb in this type 
of flow occurs near Li we expect a significant modification 
of this result because of the presence of the secondary star 
(this effect was less serious for the IP equilibrium (^ as the 
angular momentum exchange there occurs much closer to 
the white dwarf). We therefore need numerical simulations 
using the full Roche potential. 

Clearly the condition ( pT| ) favours larger values of the 
magnetic moment /ii. However this must not be so large 
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that the system synchronizes and becomes an AM Her sys- 
tem. This condition can be expressed as (e.g. King, Frank 
& Whitehurst, 1991) 



3 A CONTINUUM OF SPIN EQUILIBRIA 

The equihbrium conditions ^ and suggest that there 
exists a continuum of equilibrium states in which 



^iM2 2-KK.m 

a3 Porb 



(14) 



R ■ < R <h 



(16) 



where /i2 is the secondary dipole and M is the mass transfer 
rate. In practice this condition requires 



^ii<10^'' Gcm^ 



(15) 



Spin evolution calculations similar to those of Wynn and 
King (1995) show that the spin period of the white dwarf 
in EX Hya reaches an equilibrium which is in excellent 
agreement with the observed Pspin/-Porb = 0.68 ratio, for 

(Ll)Adyn(Ll 



Figure 5 shows this set of equilibrium states for various val- 
ues of tmag(ii)/tdyn(ii), in the casc of a binary with system 
parameters identical to those of EX Hya. Using (^ to pro- 
vide an estimate of /ii from tmag(ii) it can be seen that the 
equilibrium states between the usual IP equilibrium 

Pspin/-Porb ~ 0.1, 5 X 10^° G cm^* </ii < 10^' G cm^ (17) 



and the EX Hya equilibrium 



<l/3 (see secZn Tand figure 5). Spin ^-'"/^'o.b - 0.68, 10^^ G cm^ < < 10^^ G cm^ 



(18) 



equilibrium was attained on a time-scale 



10^ 



yrs. 



much shorter than the orbital evolution time-scale. This 
numerical result contrasts with the simple estimate above 
which predicted a ratio of Papin/Porb =0.54 from the white 
dwarf potential alone. We conclude that (^2|) indeed spec- 
ifies a spin equilibrium state, and that EX Hya is close to 
this state. 

Figure 2 shows a sequence of views of the accretion flow 
in the equilibrium spin state of our model system. This sug- 
gests that accretion of mass and angular momentum by the 
white dwarf occurs only for a restricted range of spin phases, 
essentially when the dipole axis points towards L\ . At other 
phases the condition (^ shows that matter trying to ac- 
crete along field lines will be centrifugally ejected. Much of 
this ejecta will be re-accreted by the secondary, transferring 
angular momentum from the white dwarf back to the binary 
orbit. In equilibrium these angular momentum flows balance 
over each orbit. We note here that recent spectroscopic ob- 
servations (Wynn, Wheatley and Maxted, 1999) show strong 
evidence for just this sort of flow pattern in EX Hya. Figure 
3 shows the fractions of mass transferred instantaneously 
through Li which are accreted, returned to the companion 
star, or ejected entirely from the system for the simulation 
depicted in figure 2. It can clearly be seen that most of the 
transferred mass is accreted by the white dwarf, and this 
accretion occurs in two 'bursts' each beat cycle. It should 
be noted however, that the simulation takes no account of 
the modulation of the accretion flux by the white dwarf spin 
which may well dominate the observed light curve. 

The spin equilibrium described here is quite robust and, 
in particular, may be achieved with a large spread in blob 
parameters: only ^ 25% of the accreting mass need be mag- 
netically dominated (i.e. condition |ll] holds) to achieve the 
EX Hya equilibrium. Figure 4 shows the results of a sim- 
ulation with a greater spread in blob parameters. In this 
figure it is clear that the accretion flow consists of 2 compo- 
nents, one of which is similar to that described above and 
the other is a azimuthally symmetric distribution of parti- 
cles surrounding the white dwarf. This complex system could 
mimic many of the observational properties expected of ac- 
cretion via a disc and via a stream, whilst still attaining an 
equilibrium Pspin/Porb ratio very close to that observed for 
EX Hya. 



(with the exact value of each depending upon the mass ratio 
of the system) are defined by 

0.1 < Pspin/Porb £0.68, 10^' Gcm^</ii<10^^ G cm^ (19) 

For /ii < 5 X 10^*^ G cm^ a standard Keplerian accretion disc 
would be expected to form, and for ^ 10^* G cm^ the 
system would synchronize and be observable as an AM Her 
star. Primary magnetic moments ;$ 10"^^ G cm^ are much 
lower than those normally inferred for IPs above the period 
gap (>10^^ G cm''). Therefore, we would expect IPs above 
the period gap to evolve into EX Hya-like systems or AM 
Hers. It should be noted that each system would have a 
unique curve similar to that presented in figure 5, which 
would be a function of Porb, MilM\ and M. 

The spin period of the white dwarf in the short period 
magnetic CV RX1238-38 (Pspin/Porb — 0.4) may represent 
one of these continuum equilibrium states. Confirmation of 
this is dependent on the (as yet undetermined) mass ratio 
of the system. However, the Pspin/Porb ratio is clearly too 
high to represent the usual IP equilibrium ^. If Pspin was 
determined by Pco ~ h then (|l^) gives a rough estimate of 
the mass ratio of the system as M^jMx ~ 0.65. This mass 
ratio is rather high for a system below the period gap and 
the most likely possibility is that RX1238-38 represents one 
of the possible spin equilibria defined by (fl6|). 



4 DISCUSSION 

The reasoning above shows that EX Hya represents a new 
equilibrium spin state for magnetic CVs characterized by the 
approximate equality (|l2|) . We can now show that the nec- 
essary conditions w\ [15^ generally confine this equilibrium 
state to short orbital periods, in precise agreement with ob- 
servation. In Figure 6 we plot tmag(il)/tdyn(il) versus Porb 

for three values of (We have assumed secular mean values 
for M2,M, cf e.g. King, 1988.) Evidently the required con- 
ditions for the EX Hya equilibrium hold only for a narrow 
range 

2 X 10^^ Gcm^ < Ati < 10^" Gcm^ (20) 

(in agreement with Figure 6), and then only for 

Porb < 2 hr, (21) 

i.e. periods below the gap. For ( |ll] ) to hold at periods above 
the gap would require values of probably violating (|l5|). 
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i.e. such systems are generally synchronous (AM Her sys- 
tems). Thus we expect most systems in this equilibrium state 
to be below the period gap. 

There is a continuum of equilibrium spin states which 
connect the usual IP equilibrium and that represented by 
EX Hya. The spin equilibrium for a given Porb, M2/M1 and 
M is determined by the strength of the primary's magnetic 
moment. The other short period magnetic CV RX1238-38 
most likely occupies one of these states. The relatively low 
values of /ii ( > 10'^^ g cm^) required to attain such equilibria 
at short orbital periods ( < 2 hr) indicates that magnetic CVs 
above the orbital period gap will evolve toward long spin 
periods (either fspin > O.lPorb, or Pspin = Porb) below the 
gap. It is possible that some of these longer spin equilibria 
may be attainable above the period gap. 

We may now crudely classify magnetic CVs according 
to Porb and /Ui (Fig. 9). We see that EX Hya-type systems 
have field strengths similar to typical IPs above the period 
gap and comparable with the weakest-field AM Her systems 
below it. Systems with larger /ii are evidently AM Her sys- 
tems below the gap even if not above it. The Figure prompts 
the obvious question as to why there are so few asynchronous 
magnetic CVs below the period gap. There are a number of 
possibilities: 

(i) the field distribution has a lower bound, and such sys- 
tems are rare. DQ Her would be an example above the period 
gap. 

(ii) Such systems do exist, but are difficult to detect, par- 
ticularly below the period gap. 

(iii) Orbital angular momentum losses associated with 
the ejected material drive the orbital evolution of short pe- 
riod systems significantly more quickly than gravitational 
radiation alone. 

In the second case, we should remember that IPs are usu- 
ally found through their medium-energy ( 2 kcV) emis- 
sion. The lower mass transfer rates below the period gap 
may simply make this too weak for detection. 

We have shown that the spin period of EX Hya cor- 
responds to a new equilibrium state in which the corota- 
tion radius is comparable to the distance to the inner La- 
grange point. Moreover this equilibrium is only likely below 
the period gap, because the magnetic fields required above 
the gap would make the systems synchronous. The accretion 
flow (Fig. 1) characterizing this state has distinctive features 
which should allow observational tests. We note that Wynn, 
Wheatley and Maxted (1999) present data which provide 
strong evidence for exactly this form of accretion/ejection 
flow in EX Hya. 
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Figure Captions 

Figure 1 Numerical simulations of the accretion flow in 
EX Hya, adopting the magnetic drag prescription of King 

(1993), and a dipolar magnetic field geometry. The upper 
panels depict the system as viewed from above the orbital 
(x-y) plane, where the motion of the secondary is in a anti- 
clockwise sense and the z-axis is orientated out of the plane 
of the paper. The lower panels depict the system as viewed 
from the x-z plane. In each case a half-Gaussian distribution 
for k was used: to reflect a stream of blobs with a dense core 
and tenuous wings (higher k values). The distribution mean 
and standard deviation were tmag(I/i)/fdyn(iji) = 0.25 and 
0.025 (tailing to lower values) respectively. The left-hand 
panels show the results of a simulation which also includes 
an equal number of particles with very low k values, in or- 
der to demonstrate the robust nature of the spin equilibrium 
(see section 2). 

Figure 2 A sequence of views of the accretion flow in the 
equilibrium state for a simulation with magnetic parameters 
equal to those of the right-hand panels of figure 1. The left 
hand panels in the figure show the x-y projection of the sys- 
tem, while the right-hand panels show the corresponding x-z 
projection. Each set of panels are separated in time by 0.2 

Porb- 

Figure 3 Fractions of mass transferred through L\ which 

arc accreted by the white dwarf (solid curve) , recaptured by 
the secondary star (dot-dash curve), and escaping from the 
binary (dashed curve), for the simulation shown in figure 2. 

Figure 4 As figure 2, but for a simulation with magnetic 
parameters equal to those of the simulation shown in the 
left-hand panels of figure 1. 

Figure 5 Equilibrium Pspin/Porb values as a function of 
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tmagiLi) /tdyn{Li) , and estimated /ii for EX Hya. 

Figure 6 imag(ii)/tdyn(-£'i) verses Porb- From top to bot- 
tom the panels plot the time-scales for fj,i = 10^^, 2 x 10^^, 
4 X 10*^ G cm^. 

Figure 7 Schematic classification of magnetic CVs accord- 
ing to fii and Porb- Here DQ Her systems are taken to rep- 
resent systems which accrete via a Keplerian disc and have 

Pspin/Porb « 0.07. 
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This figure "figurel.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9909312vl 



This figure "figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9909312vl 
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This figure "figure4.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9909312vl 
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